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Abstract 



The Generalized Lotka-Volterra (GLV) model: 

Wi(t + 1) = Xwi(t) + aw(t) — cw(t)wi(t) , i = 1, ,N 

provides a general method to simulate, analyze and understand a wide 
class of phenomena that are characterized by power-law probability distribu- 
tions: 

P(w)dw ~ vj-^dw (a > 1) 
and truncated Levy flights fluctuations L a (w). 

I. INTRODUCTION 

Many natural and man-made phenomena are known to involve power-law probability 
distributions (e.g. Pareto 1897; Zipf 1949; Mandelbrot 1961, 1951, 1963; Atkinson and 
Harrison 1978; Bak et al. 1997; Bouchaud et al. 97; Cahalan and Joseph 1989; Mantegna 
and Stanley 1994, 1995, 1996, 1997; Stanley et al. 1995; Sornette et.al. 1997; Zanette and 
Manrubia 1997, Zhang et al 1997). 



*To appear in Econophysics 97, (Kluver 1998), eds. Imre Kondor and Janos Kertes 
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Power-laws are common in systems composed of units that have no characteristic size, 
and in systems made of auto-catalytic elements (Yule 1924; Champernowne 1953; Simon 
and Bonini 1958; Ijiri and Simon 1977, Anderson 1995). Moreover it has been shown that in 
systems which are not separable into energetically independent parts, the power laws take 
naturally the place of the usual exponential (Boltzmann) distribution (Tsallis 1988). 

It has been shown both theoretically (Solomon 1998; Solomon and Levy 1996) and using 
numerical simulations (Biham et al. 1998) that the Generalized Lotka-Volterra (GLV) model 
produces power-law distributions, but doesn't suffer from the problems and limitations of 
previous similar models. Since the GLV requires only a few, and not too restrictive pre- 
conditions, it may be expected to be widely applicable. 

The GLV involves three scalar parameters, and a probability distribution, each having 
clear roles in the model's interpretation. Certain model properties are "universal" i.e. inde- 
pendent of some of the model's parameters. This conceptual simplicity makes the GLV easily 
adaptable to different systems, and makes it possible to use it as a descriptive, understanding 
and comparison tool. 

The features listed above make the GLV a candidate for a general method to simulate, 
analyze and understand a wide class of phenomena which are characterized by power-law 
probability distributions and multiscale fluctuations. 

The aim of this paper is to provide a practical introduction to GLV modeling. We describe 
the basic theory, how to construct a numerical simulation of the GLV and applications to 
practical interesting systems. 

For didactic reasons we will present a few previous and simpler models before introducing 
the GLV. This will nicely partition the topics involved into more manageable parts. The 
models described in this paper are: 

• Single-agent random multiplicative process without barrier eq. (|7|). 

• Single-agent multiplicative process with fixed lower barrier eq. fl20|) . 

• Multiple agent process with barrier coupled to the mean eq. (f|0]-[|^). 
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• Single-agent linear (stochastic multiplicative/additive) process eq. fl50D 

• Generalized Lotka-Volterra (autocatalytic competing) agents eq. (|52|). 

II. SOME BASIC CONCEPTS 

Microscopic representation or agent-oriented simulation is a relatively new way of gen- 
erating and expressing knowledge about complex systems (see Solomon 1995 for a physicist 
oriented review and Kim and Markowitz 1989 for an early pioneering work in finance and 
ZhangWB 1991 for a synergetically inspired view). 

One idea this document tries to convey is the basic difference between standard methods 
of explanation, which are based on a global parameterization of a macroscopic dynamics and 
the "microscopic representation" or "agent-oriented" simulation method of explanation. 

To illustrate the older methods, think of the high-school textbook problem of computing 
the trajectory of a stone under the influence of gravitation. The stone is treated as a 
point mass (that is the global parameterization), and we use Newton's Laws (that's the 
macroscopic dynamics). 

This kind of explanation is not always useful. For instance in order to study the cracking 
and crumbling of the stone under pressure, one would have to consider its structure made of 
clusters of (of clusters ...) of smaller stone parts. In many physical, economic etc situations, 
we find it is most natural to view the system under study as a collection of "microscopic" 
similar units (or "agents"), which interact among themselves in some well-defined way and 
(co-)evolve in time. 

Experience have shown that even a simple dynamics applied to a system consisting of 
many similar interacting agents may produce interesting macroscopic effects. The proper- 
ties which arise from the collective behavior of many similar elements are called emergent 
properties. 

Some examples of systems with emergent properties are: 
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Animal (and human) populations (cities/countries) are composed of individuals which 
are constantly born (and die), and compete with other members of their species (for food, 
mates etc). 

The classical (scalar) Lotka-Volterra system: 



was invented specifically in order to model the time evolution of such populations w(t). 
We will show that its straightforward multi-agent generalization eq. (2) which treats the 
population as a collection of sub-populations Wi, presents emergent properties very different 
from the eq. (1) which treats the entire population as a single variable (similarly to the 
stone/point mass example above). 

The stock market is parametrized macroscopically by the index w(t) which is proportional 
to the capitalization (total worth of the money invested) in the traded equities. For a fixed 
amount of shares, this is a measure of the price of the shares traded in the market. Models 
similar to the stone/point-mass/Newton law example above were introduced in the past 
(Baillie and Bollerslev 1990) and lead to models which treat the index/share price w(t + 1) 
as a single entity with time evolution governed by stochastic differential equations. 

The microscopic representation method models the market index w(t) as an collective 
quantity emerging from the interactions of a macroscopic number of traders. 

More precisely, the "microscopic representation" of the stock market is composed of 
many investors i = 1, N each having a certain personal wealth Wi(t). The investors buy 
and sell stocks, using more or less complicated strategies designed to maximize their wealth. 
The resulting index is then proportional to the total sum of the individual investors wealth. 

Another ("orthogonal" ) way to characterize the financial markets is to consider as elemen- 
tary degrees of freedom the individual stocks (i.e. the capitalization Wi(t) of the individual 
companies % traded in the market). 

We are going to use the Generalized Lotka-Volterra (GLV) equation system 




(1) 



Wi{t + 1) = X(t)wi(t) + a(t)w(t) - c(t)w(t)wi(t) 



,N 



(2) 
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to construct a very simple model of the stock market. In spite of being very simple, the 
model yields (and in doing so - explains) the Pareto distribution of wealth, an economical 
fact left for a long time (100 years!) without satisfactory explanation. In the case of the 
previous example of populations/cities/countries sizes, the system eq. (0) results from the 
assumption that the Wi individuals belonging to each city/country % interact similarly (in the 
stochastic sense) to all the individual members of the other cities/countries. More precisely, 
eq. @ would result from the following assumptions: 

- for each of the Wi citizen of a city, there is a A probability that he will attract (or give 
birth to) a new citizen. One may include as a negative contribution to A the probability for 
the citizen to die or leave the town without preferred destination. 

- each of the Wi citizens of a city has a probability c/N to be attracted by one of the Nw 
citizens of the existing cities and leave the town i. 

- the citizens which left own town but are not bound to a particular town, have the 
probability a to end up in any of the i towns. 

The analysis below shows that if the last term is negligible, the exponent a is likely to 
be 1. 

Systems made of auto-catalytic (terms A and a) competing (term c) elements are wide- 
spread in nature. Another example, very different from the previous: clouds are collections 
of water droplets, which grow by taking in water vapor from the surrounding air. Agent- 
oriented models may explain the observed distribution of cloud sizes and shapes without 
going into all the very complicated physical details. 

To make things more concrete, we will often use in the following discussion terminology 
borrowed from the GLV model of the stock market. However, our results and insights will 
be quite generic and applicable to other systems. 

For example, we will often use the term 'traders' i instead of microscopic elements and 
talk about the 'wealth' Wi of each trader and the "average wealth" w(t): 

™(*) = 4X>(*) (3) 

i 
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All these names should be suitably renamed when used in other contexts. 

III. PROBABILITY DISTRIBUTIONS AS PREDICTIVE OUTPUT OF 

MODELING 

Modeling the stochastic evolution in time t = 1, ,T of a system of type eq. ([!]) 

produces a lot of raw data. It is not very useful to produce more and more sequences of 
numbers 

w{t); t = l,....,T 

without a way to analyze this data. 

A natural way is to collect a lot of such sequences and analyze them statistically. A way 
to look at it is to imagine a large set of (uncoupled) traders Wi, i — 1, N, subject to the 
same stochastic dynamics, and compute their distribution. More precisely, assuming each 
data sequence i produced by eq. (P describes the evolution of one trader i 

Wi(t); t = l,....,T 

we may compute over a large set of traders i = 1, N the number of traders N t (w) with 
wealth Wi(t) in the interval (w, w + dw). In the limit of infinite N this defines the individual 
wealth probability distribution 

P t (w)dw = N t (w)/N (4) 

This way of looking at the problem is conductive to a model like eq. (0) in which the 
various traders do interact e.g. through terms including their average w. As we will see, 
this is a key ingredient in solving some of the problems with the one-agent models of the 
type eq. (|TJ) and of course a way to make the models more realistic. 
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IV. THE GENERAL SHAPE OF THE PROBABILITY DISTRIBUTIONS 

It may be a little surprising, but in certain conditions, even for non- stationary interacting 
non- linear dynamical systems of the type eq. (0), we may predict (without performing any 
computer runs) how the wealth distributions Pt(w) produced by our models will look like. 

Before going into a more detailed analysis we list here some obvious properties which 
the probability distribution has to obey 

• The distribution vanishes for negative values. 

P(w) = for w < 

This property follows from the way we choose the initial wealth values, and the na- 
ture of the dynamics. In practice this corresponds to the fact that population sizes, 
capitalization of companies, cloud sizes, etc. assume only positive values. 

• The value zero usually has vanishing probability distribution P(0) = 0. This follows 
from the continuity of the distribution, and the previous characteristic. 

• The distribution has a tail which may obey a power-law 

P(w) ~ uT 1_a (5) 

or not. In any case, distributions must decay to zero at infinity 

lim P(w) = 

W— >oo 

otherwise the probability to have finite Wi values would vanish. 

• The distribution has a maximum. A "nice" continuous and non-negative function will 
have a maximum between two zeroes. Here we have one at zero, and one at infinity. 

The remaining non-trivial issues are: 
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• The behavior of the tail, in particular whether it decreases as a power-law eq. (|5|), 
log-normal, exponential or some other shape. 

• The location and height of the maximum 

• The behavior of the rising part 

• The relation between the tail shape and the parameters of the model. We will see that 
in certain conditions, the tail behavior is quite universal in as far as it is quantitatively 
independent on most of the parameters and it is preserved even in generically non- 
stationary conditions. 

In this article we will be interested mainly in the shape of the tail of the probability 
distribution P(w), in particular whether it is a power-law eq. (|5|) or not. 

The importance of the power laws is far from purely academic: the power laws are the 
bridges between the simple microscopic elementary laws acting at the individual level and 
the complex macroscopic phenomena acting at the collective level. They insure that the 
dynamics does cover many dynamical scales rather than acting only at the smallest and/or 
largest scales of the system. 

Moreover, power laws insure that the dynamics of the intermediate scales of the system 
is largely independent on the microscopic details and the macroscopic external conditions 
constraining the system. Consequently the macroscopic complexity is not the direct and lin- 
ear result of the microscopic details but rather a generic consequence of the self-organization 
taking place in systems with many interacting parts and inter-related feedback loops. 

This means that one can hope to extract macroscopic and mesoscopic laws which hold 
for large classes of microscopic laws in an universal way. 

This unifying and predictive power of the power laws made them the central objects 
in many branches of science starting with quantum field theory and statistical mechanics 
and ending with ecology and economics. In fact the necessity for such laws in fundamental 
physics lead to very "un-natural" (t'Hooft 1979) fine-tuning procedures designed to enforce 
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them. 

One can say that the emergence of the power laws is a sine qua non condition for the 
emergence of the macroscopic world out of local microscopic elementary laws of nature. As 
such, it becomes a fundamental natural law on its own. 

The present paper shows that power laws emerge generically in the most 
simple and natural models which were considered in the past in the modeling 
of chemical, biological and social systems. 

V. MULTIPLICATIVE RANDOM DYNAMICS AND LOG-NORMAL 

DISTRIBUTIONS 

Power-law probability distributions eq. fl5|) with exponent —1 (i.e. a = 0), can be 
obtained analytically from a multiplicative process (Redner 1990; Shlesinger 1982): 

w(t+ 1) = X(t)w(t) (7) 

where the random variables A(t) are extracted from a fixed probability distribution 11(A) 
with positive support. 

Indeed, in order to obtain the distribution of w(t) in the large t limit, one takes the 
logarithm on both sides of eq. (|7|) and uses the notations ji = In A, x = Inw. With these 
notations eq. (^) becomes: 

x(t + 1) = n(t) + x(t) (8) 

The respective probability distributions p(/i) and V(x) for fi and x are related to the distri- 
butions for A and w by the identities: 

p(lnAMlnA) = U{X)d\ (9) 

and respectively 

V(\nw)d(\nw) = P(w)dw (10) 
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which mean: 



p(/jL)dfjL = (exp /x)n(exp p)dp (11) 

and respectively 

V{x)dx = (exp x)P(exp x)dx (12) 

The interpretation of eq. (||) is that x(t + 1) is the sum of the constant x(0) = lnw(O), 
and t random variables fi(t) = lnA(t) extracted from the fixed distribution p(/x). Under the 
general conditions of the Central Limit Theorem (CLT) we get for x(t) at large t the normal 
distribution: 

V t (x) ~ exp - l 2 ; (13) 

where 

< x >= tv = t < In A > (14) 

and 

a 2 iA =< (hiA) 2 > - < In A > 2 (15) 

Note that the width of the normal distribution Vt(x) (the denominator of the expression 
in the exponential in eq. (0)) is 

°l = t°lx (16) 

and increases indefinitely with time. This means that the distribution Vt(x) becomes inde- 
pendent of x 

V t {x) ~ \ (17) 

V 27r0r lnA^ 

in an ever- increasing neighborhood of < x > (where the exponential is close to 1). 
Transforming eq. fll7D back to the w variables by using eq. fllOf) one gets 
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P(w)dw = rf(lnw) (18) 



i.e. 



P(w)dw ~ efa; (19) 

Graphically, this means that as the time goes to infinity, P(w) is an ever "flattening" 
distribution approaching w~ l in an ever expanding neighborhood. 

VI. NON-INTERACTIVE MULTIPLICATIVE PROCESSES WITH FIXED 

LOWER BOUND 

In order to obtain a power-law probability distributions P(w) ~ w^ 1 ^ with exponent 
a > the multiplicative process eq. (0) has to be modified (Yule 1924, Champernowne 
1953, Simon and Bonini 1958, Ijiri and Simon 1977) in such a way that the variation of w(t) 
under eq. (0) will be constrained by a lower bound (or barrier) w m i n 

w(t) > w min . (20) 

In terms of x(t) = lnw(t) , eq. (8) becomes supplemented by the lower bound: 

x(t) > x min = ln(uw) (21) 

More specifically the dynamics of w(t) consists at each time t in the updating eq. (^) 
(or - equivalently - eq. @) except if this results in w(t + 1) < w m i n (or, equivalently 
x(t + 1) < x min = \n(w min )) in which case the updated new value is w(t + 1) = w min 
(respectively x(t + 1) = x min ). 

These modifications are obviously not allowed by the Central Limit Theorem, and conse- 
quently, the derivation in the previous section which led to the log-normal distribution eqs. 
(13) (19) cannot be applied. 

Instead, one can obtain intuition on the modified system eq. (7) (20) by realizing that in 
fact the system eq. (f|) with the constraint eq. (^1) can be interpreted as the vertical motion 
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of a molecule in the earth gravitational field above the earth surface. In this case the "earth 
surface" is x m i n , the gravitationally induced downward drift is v =< x > jt =< In A > and 
the diffusion per unit time is parametrized by the squared standard deviation af nX . 
For v < this is the barometric problem and it has the static solution: 

P(x) ~ exp-x/kT (22) 

which when re-expressed in terms of u>'s becomes (cf. eq. (10)): 

P{w)dw~e- {lnw)/kT dlnw (23) 

i.e. 

P(w) ~ w- l - llkT (24) 
In order to estimate the value of the exponent 

a = 1/kT (25) 



one can substitute the form eq. ( |24|) in the master equation which governs the evolution of 
the probability distribution of the process eq. (|7|). The master equation expresses the flow 
of probability between the various values of w as w(t) is updated to w (t + 1) = Xw(t). More 
precisely, the probability for w(t + 1) = \w(t) to equal a certain value w is the integral over 
A (weighted by 11(A)) of the probabilities that w(t) = w/\. This leads in equilibrium to the 
relation: 

P(w) = J U{X)P(w/X)d(w/X) (26) 
i.e., substituting (24) (25) into eq. (26): 

w- 1 -" = y n(A)(w/A)- 1 - Q rf(w/A) (27) 
or, by dividing both members by w^ 1 ^": 

1 = J X a U(X)dX (28) 
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This means that the value of a is given by the transcendental equation: (Solomon and Levy 



One may wonder why this equation does not hold for the problem in the previous section 
(the system eq. (11) without the lower bound eq. (20)) since the lower bound w m i n does 
not seem to appear in the formula eq. fl29P . The answer is that in the absence of the lower 
bound, there is no stationary equation (|2q) as the system evolves forever towards the non- 
normalizable solution P(w) ~ w' 1 . In fact, formally, a = is a solution of eq. (^9|). This 
"solution" is in fact the relevant one in the case < In A >> 0. 

Another nontrivial fact is that the limit w m i n — > leads to a — > 1 rather than a — > 
(which is the value in the total absence of a lower bound). This nonuniform behavior of 
the limits iV — > oo and w m i n — > is related to the non-trivial thermodynamic limit of the 
system (Biham et al. 1998). 

One of the problems with the solution to eq. ( |2"9"D is that while it is independent on w min ,it 
is highly dependent on the shape and position of the distribution 11(A) of the random factor 
X(t) and consequently it is highly dependent on the changes in the dynamics eq. (11). 

In the systems introduced in the following sections, the situation will be the opposite: 
the characteristics of 11(A) will be largely irrelevant and the lower bound eq. (20) will play 
a central role. 

In particular, the distribution of the social wealth P(w) and the inflation (dw/dt) of the 
system will depend on the social security policy, i.e. on the poverty bound w min below of 
which the individuals are subsidized (Anderson 1995). The relevant parameter is the ratio 



1996): 



< \ a >= 1 



(29) 



q = w min /w 



(30) 



between the minimal wealth w. 



rmn 



and the average wealth: 



w = l/N y £w i 



(31) 
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In the presence of inflation, it is q rather then w m i n which has to be fixed since a fixed minimal 
wealth independent on the changes in the average wealth would be not very effective for long 
time periods. 

To prepare the formalism for those more sophisticated applications we will deduce below 
an identity relating the exponent a = 1/kT to the ratio q eq. (30). 

Such a relation can be obtained using in eq. (p4|)(25) the fact that the total probability 
is 1: 

/oo 
P{w)dw = 1 (32) 

i.e. 

Const. / w~ l - a dw = 1 (33) 
and the fact that the average of w is w): 

wP{w)dw = w (34) 
i.e 
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Const. / w a dw = w (35) 
By extracting Const from eq. (j33[): 

Const = aw^ in (36) 
and introducing this value in eq. (|35|) one gets the relation: 

Kml-w™/! 1 -a)]=w (37) 
By dividing both members by w-^^ one obtains: 

- w min /w = l/a-l (38) 



Considering the definition eq. (|30D this reads: 
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« = l/(l-g) (39) 

This is a quite promising result because for the natural range of the lower bound ratio 
< g < 1/2 it predicts 1 < a < 2 which is the range of exponents observed in nature (which 
is far from the log- normal value a = 0). 

The main problem with the power-low generating mechanism based on the single-agent 
dynamics eq. (11) (20) is that it works only for negative values of < In A > (the barometric 
equation does not hold for a gravitational field directed upwards). 

This means that (except for the neighborhood of w m i n ) w(t + 1) is typically smaller than 
w(t). This is not the case in nature where populations, economies are expanding (at least 
for certain time periods). 

Moreover, the exponent a of the power law is highly unstable to fluctuations in the 
parameters of the system. In particular, trying to model large changes in w one is lead to 
large fluctuations of q and consequently (cf. eq. (|3~9"D) to large variations in the exponent 
a. This again is in disagreement with the extreme stability of the exponents a observed in 
nature. 

We will see later how the multi-agent GLV solves these problems. The key feature is to 
introduce interactions between the individual elements Wi through the inclusion of terms (or 
lower bounds) proportional to w. 



VII. MULTIPLICATIVE PROCESSES COUPLED THROUGH THE LOWER 

BOUND 



In the previous section we showed that power-laws eq ([24]) can be obtained from multi- 
plicative stochastic dynamics with lower bound in the same way that the exponential laws 
eq. (22) can be obtained in additive stochastic dynamics eq. (8) bounded from below eq. 
(21). 

The problematic points were: the instability of the exponent a to variations in the 
average wealth/population and the fact that the mechanism accounts only for "deflat- 
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ing" / "shrinking" of w. 

However, the solution of these difficulties is already apparent in eq. we consider 

a system of iV degrees of freedom wi (with % — 1, ...N) which is governed by the following 
dynamics: 

at each time t, one of the iWj's is chosen randomly to be updated according to the formula: 



while all the other w^s are left unchanged. At each instance, the random factor A is extracted 
anew from an i-independent probability distribution n(A). 

The only exception to the prescription eq. (fHf) is if following the updating eq. (|40D 
Wi(t + 1) (or other w^s) end up less then a certain fixed fraction < q < 1 of the average 



Wi(t + 1) = \(t) Wi (t) 



(40) 



Wi(t + 1) < qw{t) 



(41) 



The prescription, if eq. (pEID happens, is to further update the affected lOj's to: 



Wj(t + 1) = qw(t). 



(42) 



The wonderful property of the system eq. (f40|-[42|) is that when re-expressed in terms of the 



variables 



Vi(t) 



Wi(t)/w(t) 



(43) 



it leads to a system very close to the system eq. Q7D , (fZO]) : 



Vi(t+1) = X(t)vi(t) 



(44) 



Vi(t + 1) > q 



(45) 



where the effective multiplicative factor in eq. ([0]) is: 



\{t) = \(t)w(t)/w(t+l) 



(46) 
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Note that this solves the problem of < In A >> because the multiplication with the 
( "renormalization" ) factor w(t)/w(t + 1) takes care that the v distribution is never running 
to infinity (in fact it insures that v(t) = 1 always). 

Another way to see how this solves the problems with the single "particle" dynamics 
eqs. flTDflSOP is to realize that now, with a qw lower bound, even if the average w(t) runs to 
infinity, the lower bound runs after it in such a way as to insure a stationary value of a. 



Indeed, according to eq. (|24[) , the Vi dynamics eq. (fHHiq) leads to a V{ distribution 

P(v) ~ v- l - a (47) 

which in turn implies 

P(w) ~ w- 1 -* (48) 

Since we do not have a close analytic formula for A eq. (46) in terms of the model parameters, 
the transcendental equation eq. (p9|) is not useful here. 

However, since v = 1 by definition, one gets according eq. ( |39"D the following close formula 
for the exponent a: 

a = 1/(1 - q) (49) 

This means that even for significantly time- varying distributions lit (A) in eq. (40), the 
exponent of the power law remains time-invariant. 

The above results are quite non-trivial in as far as they predict the stochastic 
behavior of highly interactive and time non-stationary systems eqs. (40)-(42) by 
relating them formally to non- interacting static statistical systems eqs. (7) (20). 

We will continue this line in the following sections an reduce the highly nontrivial GLV 
system to a rather simple single-agent linear stochastic equation. 

VIII. THE SINGLE- AGENT LINEAR STOCHASTIC EQUATION 



The discrete equation below is a more elaborate version of the single-agent model intro- 
duced in Section 6. Here the lower bound is supplied effectively by an additive random term 
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p, instead of an explicit barrier eq. (20): 

w(t + 1) = X(t)w(t) +p(t) (50) 

The equation ( |50"D may (crudely) describe the time-evolution of wealth for one trader in the 
stock market. In this context w may represent the wealth of one trader, who performs at 
each time-step a stock transaction. w(t) is the wealth at time t, and w(t + 1) is the wealth 
at time t + 1. 

A and p are random variables extracted from positive probability distributions. A may 
be called the "success factor", and p is a "restocking term", which takes into account the 
wealth acquired from external sources (e.g. using the state built infra- structure, subsidies 
etc). 

The dynamics eq. (|50|) leads to a power-law distribution P{w) ~ w~ 1 ~ a , if < In A >< 0. 
More precisely, the P(w) has a power-tail for values of w for which p(t) is negligible with 
respect to X(t)w(t). If < In A >> the distribution is a log-normal expanding in time, which 
in the infinite time limit corresponds to a power-law with exponent — 1 — a = — 1. I.e. the 
mechanism eq. (50) does not explain expanding economies with a > 0. 

Since the distributions (and initial values of w) are positive, the p(t) term keeps the 
value of w(t) above a certain minimal value of order p. Therefore, for large enough values 
of w the dynamics is indistinguishable from our previous model, the multiplicative process 
(7) with fixed barrier eq. (f2"0p. 

It is therefore not surprising that it can be rigorously proven (Kesten 1973) that eq. (|50|) 
leads to a power law eq. fl24|) with exponent a given by the transcendental equation (|29|): 

< X a >= 1 (51) 

It is again notable that (as in the case of the w m i n lower bound), the exponent a is totally 
independent on the distribution of the additive term p while it is highly sensitive to the shape 
and position of the 11(A) distribution. We will see that in GLV the situation is reversed. 
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IX. THE GENERALIZED LOTKA-VOLTERRA SYSTEM 



The GLV (Solomon and Levy 96) implements a more complex dynamics than the eq. 
(|50|). The advantages however overweight the difficulties: 

• The GLV model insures a stable exponent a of the power-law P(w) ~ w~ l ~ a even in 
the presence of large fluctuations of the parameters. 

• The value of < A > and the average w can vary during the run (and between the runs) 
without affecting the exponent a of the power-law distribution. 

• A can take typical values both larger or smaller than 1. 

The GLV is an interactive multi-agent model. We have N traders, each having wealth 
Wi, and each Wi is evolving in time according to: 

Wiit + 1) = X(t)wi(t) + a(t)w{t) - c(t)w(t)wi(t) (52) 

Here w is the average wealth, which supplies the coupling between the traders: 

w = (wi + w 2 + . . . + w N )/N (53) 

A is a positive random variable with a probability distribution 11(A) similar to the A success 
factor we used in eq. (50) in the previous section. The dramatic difference is that now, A 
can take systematically values larger than 1 and in fact its distribution can vary in time 
(and have time intervals with < In A > both smaller and larger than 0). 

The coefficients a and c are in general functions of time, reflecting the changing conditions 
in the environment. 

The coefficient a expresses the auto-catalytic property of wealth at the social level, i.e. it 
represents the wealth the individuals receive as members of the society in subsidies, services 
and social benefits. That is the reason it is proportional to the average wealth. 

The coefficient c originates in the competition between each individual and the rest of 
society. It has the effect of limiting the growth of w to values sustainable for the current 
conditions and resources. 
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X. REDUCING GLV TO A SET OF INDEPENDENT EQUATIONS (50) 



Summing the GLV eq. (|52|) over i, and taking the local time average one gets for the 
local time average < w > (t) an equation similar to the scalar LV equation eq. (1). 

N < w >=< \ >N<w> + < a(t) > N < w > — < c(t) > N < w > 2 (54) 

which gives : 

-1+ < A > + < a(t) > . . 

< w >= — — (55) 

< c(t) > y ' 

Neglecting the fluctuations of the average < w > (t) during the updating of a single 
individual i, one can substitute < w > (t) for w it in the last term of eq. (|52"D and regroup 
the terms linear in wf. 

Wi {t + 1) = [1 + \{t)- < X(t) > - < a{t) >]wi(t) + a(t)w(t) (56) 

Introducing like we did in eq. fl4~3| ) wealth values normalized by the average wealth: 

Vi(t) = Wi{t)/w{t) (57) 

the equation (|56|) becomes: 

Viit + 1) = [1 + X(t)- < X(t) > - < a(t) >)]vi(t) + a(t) (58) 

where we neglected again the fluctuation of w during the time t— >t + 1 and put = 1. 

This can be justified rigorously for a > 1 in the large N limit because then the size of the 
largest Wi is (cf. Solomon 1998) of order 0(wN l ~ a ) « 1 and therefore the changes in w(t) 
induced by any Wi updating are negligible of the leading order. In finite systems N < oo 
and for a < 1, there are (computable) corrections. 

The key observation now is that the system eq. (^) has the form of iV decoupled 
equations of the form ( |50D with the effective multiplicative stochastic factor 

A = 1 + A(t)— < A(i) > — < a(t) > (59) 
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This means that the results of the single linear stochastic agent model eq. ([50]) can be 
applied now to the t>i(£)'s in order to obtain 

P(v) ~ v- l - a (60) 

which in turn implies 

P{w)^W- 1 ~ a (61) 

with the exponent a dictated by the equation eq. (^TT) (59) : 

< [1 + X(t)- < X(t) > -a(t)] a >= 1 (62) 

This shows that the wealth distribution in a economic model based on GLV is a power 
law. In addition, it is easy to see using (58) that the elementary steps in the time variations of 
the average wealth are distributed by a (truncated) power law. Accordingly, it was predicted 
that the market returns will be distributed by L a (w) a truncated Levy distribution of index 
a (Solomon 1998). This turned out to be in accordance with the actual experimental data 
(Mantegna and Stanley 1996). 

Note that similarly to the passage from the single agent system eq. (7) (20) to the many 
agents system coupled by the lower bound eq. (40) (41), here too, the formal reduction of 
the GLV to a single agent system implies very different properties in the actual "physical" 
system. 

In particular while in the single agent system the average of A was crucial in the fixing of 
the exponent a and in fact A had to have an average less then 1, in the multi agent model, 
the average of A cancels in the expression ( [59]) for A. On the other hand, while the additive 
term in eq fl50| ) had no role in the fixing of a, the corresponding term a is one of the crucial 
factors determining a in GLV (cf. eq. (|62])). Moreover, while the parameters in the eq. ( |50| ) 
model allowed time variations in the w only at the price of variations in the exponent a, 
in the GLV system, one can arbitrarily change the ecological/economic conditions c so as 
to vary the total population/wealth by orders of magnitude (cf. eq. (55)) without affecting 
the power law and its exponent (A is independent on c and so is the solution of < A a >= 1. 
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XI. THE FINANCIAL INTERPRETATION OF GLV 



In this section we discuss the various terms appearing in the equations, their interpreta- 
tion in the financial markets applications, their effects and their implications for the financial 
markets phenomenology. 

We first discuss the assumption that the individual investments/gains/losses are propor- 
tional to the individual wealth: 

w(t + l) = Xw(t) 

This is actually not true for the low income/wealth individuals whom incomes do not orig- 
inate in the stock market. In fact the additive term aw tries to account for the departures 
related to additional amounts originating in subsidies, salaries and other fixed incomes. 
However, for the range of wealth where one expects power laws to hold (w > w) it is well 
documented that the investment policies, the investment decisions and the measured yearly 
income are in fact proportional to the wealth itself. 

The statistic uniformity of the relative gains and losses of the market participants is 
a weak form to express the fairness of the market and the lack arbitrage opportunities 
(opportunities to obtain systematically higher gains A — 1 than the market average without 
assuming higher risks): for instance, if the distribution of A would be systematically larger 
for small-w-investors, then the large-w-investors would only have to split their wealth in 
independently managed parts to mimic that \ow-w superior performance. This would lead 
to an equalization of their A to the A of the low-w investors. Therefore in the end, the 
distribution 11(A) will end up w-independent as we assumed it to be in GLV from the 
beginning. 

Almost every realistic microscopic market model we have studied in the past shares this 
characteristics of w-independent 11(A) distribution. 

Turning to the terms relevant for the lower-bound w region, the assumption that the 
average wealth contributes to the individual wealth aw and the alternative mechanism as- 
suming a lower bound proportional to the average qw are both simplified mechanisms to 
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prevent the wealth to decrease indefinitely. 

In practice they might be objectionable: it is not very clear that the state subsidies can 
be invoked to save bankrupt investors (though this happens often when large important 
employers are at risk or when their collapse would endanger the stability of the entire 
system). In any case the term aw is appropriately taking into account the arbitrariness 
of the money denominations. More precisely, assume that (by inflation or by currency 
renaming) the nominal value of all the money in the economy becomes 10 times larger. 
Then the subsidies term aw will become 10 times larger preserving in this way the actual 
absolute value. So to speak, multiplying a quantity by w expresses it in "absolute currency". 

These mechanisms controlling the lower bound behavior of the system may be considered 
as just parametrizations of the continuum flow of investors/capital to and from the large-w 
investor ranges relevant for the financial markets trading. It is however possible that the 
subsidies to the very poorest tumble their way (through the multiplicative random walk) into 
the middle classes and end-up in the large-w scaling region in the way our models suggest. 

More theoretical and experimental research is necessary in order to discriminate between 
the various alternatives (or in order to recognize them as facets of the same phenomenon). 

For instance one can look at the relation between q and a as basically kinematic in the 
sense that given the power law, it is unavoidable that the lower bound qw would govern the 
exponent of the power law. 

It would still be interesting to study in detail how the wealth pumped at the lower-bound 
barrier makes its way to the large w tail and or, alternatively in the case of stationary w 
to find the way in which the additive subsidies to the low-w individuals in our models are 
covered by the " middle classes" . 

For instance in the qw lower-bound model, the dominant mechanism of extracting wealth 
from the middle class seems to be the accelerated inflation: the updating of the individual 
WiS induces changes in w which in turn induces changes in the position of the lower bound 
qw which leads to the necessity to subsidize immediately all the w^s situated between the 
old and the new poverty line. This in turn leads to increase in the w and to the completion 
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of the positive feedback loop. When looking at the normalized wealth vi = Wi/w, this 
accelerated inflation is effectively a proportional tax which reduces the A gains to smaller 
relative gains A = \w(t)/w(t + 1). 

In the case of the GLV model, there is no inflation (for constant c, w is - modulo local 
fluctuations - constant in time). The extraction from the middle class of the money for 
subsidies is quite explicit from the way the subsidies term aw is affecting negatively (—a) 
the A gain value eq. (39). However most of the loss in A eq. (39) is the term (— < A >) 
and is due mainly to the competition between the large traders —cWiW. This term is in 
fact equivalent with enforcing a proportional tax. The best way to see it is to recall from 
above analysis that multiplying the wealth-proportional quantity cwi by w expresses it in 
"absolute currency". 

It would appear that the government has the choice of enforcing such a proportional tax 
and keep w roughly constant, or just print the money which it dispenses to the poor and 
let the inflation tax the other classes. One can of course make compromises between these 
2 extremes by allowing both taxation and inflation. In any case, the net result can be only 
variations in the relative size of the middle class (variations of a) as the generic emergence 
of the power law will be very difficult to avoid. 

A government may try to get a more equalitarian distributions (increasing a) by in- 
creasing q. However this would bring an even larger population in the neighborhood of the 
poverty line. This would require more and more frequent subsidies to enforce and conse- 
quently (according to the analysis above) larger taxes and/or faster inflation. 

Fearing this, another government might decide for low values of a. This would not 
only mean an increase of the ratio between the richest and poorest which might be morally 
questionable but also lead to dramatically unstable fluctuations in the market (e.g. in w). 
Indeed, one can show that typically the largest trader owns 0(N 1 ~ a ) of the total wealth. For 
finite N and very low q, a may be shown to drop below 1 (contrary to eq. (49)). This would 
imply that almost all the wealth is owned by the largest lUj. It is well know however that 
the single agent discrete logistic map eq. (1) leads generically to chaotic unstable dynamical 
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regime (May 1976). 

In between these extremes it might be that there is no much choice of dynamically 
consistent a values outside the experimentally measured range 1.4 < a < 1.7. These values 
correspond (through eq. (49)) to poverty lines between 0.3 < q < 0.4. Larger values of 
q would imply that almost everybody is subsidized while a significantly lower value would 
mean the poorest cannot literally live: the average wages in an economy are automatically 
tuned as to insure that a family of 3 can fulfill the needs considered basic by the society on 
a one-salary income. Somebody earning less than will therefore have serious problems 
to live a normal life. In fact, such a person might be exposed to hunger (even in very rich 
economies) since even the prices of the basic food are tuned to the level of affordability of 
the average family in the given economy. 

Far from implying a fatalistic attitude of the economic facts of life, the analysis of our 
simplified models might lead in more realistic instances to ideas and prescriptions for dy- 
namically steering the social economic policies to optimal parameters both from the techni- 
cal/efficiency and from the human/moral point of view. 

They might help transcend (through unification) the traditional dichotomy by which 
science is only in charge with deciding true from false while humanities are in charge only 
with discerning good from bad. 

XII. FURTHER ECONOMICS APPLICATIONS OF GLV 

As mentioned a few times, the GLV systems can be applied to many power laws. 

Even within the financial framework there are a few apparently different realizations of 
the GLV dynamics and of the power laws. 

For instance, one can consider the market as a set of companies i = 1, ,N whose 

shares are traded and whose prices vary in time accordingly. 

One can interpret then Wi as the capitalization of the company i i.e. the total wealth of 
all the market shares of the company. 
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The time evolution of Wi can still be represented by eq. (0). In this case, A represents 
the fluctuations in the market worth of the company. For a given total number of shares A 
is measured by the change in the individual share price. These changes take place during 
individual transactions and are typically fractions of the nominal share price (measured in 
percents or in points). 

With such an interpretation, a represents the correlation between the worth of each 
company and the market index. This correlation is similar for entire classes of shares and 
differences in the a's of various economic sectors result in the modern portfolio theory in 
"risk premia" which affect their effective returns A in a way similar to our formulae in a 
previous section (low correlation a corresponds to larger effective incomes). 

c represents the competition between the companies for the finite amount of money in 
the market (and express also the limits in their own absolute worth). We do not need to 
consider c a constant. Time increases in the resources may lead to lower values of c which in 
turn lead to increases in w. However, as seen before, such changes do not affect the exponent 
of the power law distribution. 

In this interpretation, the GLV model will predict the emergence of a power-law in the 
probability distribution of company sizes (capitalization). In particular, this would imply 
that the weights of the various companies composing the S&P 500 are distributed by a 
power law. In turn this would imply that the S&P fluctuations follow a truncated Levy 
distribution of corresponding index. 

Yet another interpretation is to consider W{ as the size of coordinated trader sets (i.e. 
the number of traders adopting a similar investment policy) and assuming that the sizes of 
these sets vary self-catalytically according to the random factor A while the a term represents 
the diffusion of traders between the sets. Such an autocatalytic dynamics for the trading 
schools is not surprising as their decision processes involve elements similar to the use of 
common language, common values, which are central in the dynamics of languages and 
nations (which fulfill power laws) too. The nonlinear term c represents then the competition 
between these investing schools for individual traders membership. In fact if a and c fulfill 
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the relation a/c = w, the corresponding two terms in GLV taken together represent (in 
"absolute currency" w) the act of each of the schools loosing (shedding) a certain c fraction 
of its followers which are then spread uniformly between the schools. 

Such an interpretation of the GLV would predict a power-law in the distribution of trader 
schools sizes (similar to countries sizes). Again, such a distribution would account for the 
truncated Levy distribution of market fluctuations (with index equal to the exponent of the 
power-law governing trader sets sizes). 

It would be interesting to discriminate between the validity of the various GLV in- 
terpretations by studying experimentally the market fluctuations and the distributions of 
respectively individual wealth, companies capitalization and correlated investors sets. It is 
not ruled out that some of these interpretations may be consistent one with the other. This 
would be consistent in turn with the modern portfolio theory claim on the existence of a 
"market portfolio" (stochastically) common to most of the traders. 

Further developments of the GLV model will include introducing variable number of 
agents N, studying the role of the discretization of the Wi changes (due to indivisible units 
like people, shares, etc), taking into account the influence of the history of w(t) on A (like 
in the case of investing strategies, crashes memory, etc.). 
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